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Abstract. I describe the construction of a large class of asymptotically fiat initial 
data with non-vanishing mass and angular momentum for which the metric and the 
extrinsic curvature have asymptotic expansions at space-like infinity in terms of powers 
of a radial coordinate. I emphasize the motivations and the main ideas behind the 
proofs. 
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1 Introduction 

Suppose we want to describe an isolated self-gravitating system. For example 
a star, a binary system, a black hole or colliding black holes. Typically these 
astrophysical systems are located far away from the Earth, so that we can re- 
ceive from them only electromagnetic and gravitational radiation. How is this 
radiation? For example one can ask how much energy is radiated, or which are 
the typical frequencies for some systems. This is the general problem we want 
to study. These systems are expected to be described by the Einstein field equa- 
tions. It is in principle possible to measure this radiation and compare the results 
with the predictions of the equations. 

There are no explicit solution to the Einstein equations that can describe 
such systems. Since the equations are very complicated it is hard to believe that 
such explicit solution can ever be founded. Instead of trying to solve the complete 
equations at once, we use the so called 3+1 decomposition. We split the equations 
into "constraint" and "evolution". First we give appropriate initial data: that is, 
a solution of the constraint equations. Once we have chosen the initial data, the 
problem is completely fixed. Then, we use the evolution equations to calculate 
the whole space-time. From the evolution we can compute physical relevant 
quantities like the wave form of the emitted gravitational wave. This method 
of solving the equations is consistent with the idea that in physic we want to 
make predictions, that is: knowing the system at a given time we want to predict 
its behavior in the future. It is of course in general very hard to compute the 
evolution of the data, one has to use numerical techniques. The question we want 
to study here is: what are the appropriate initial data for an isolated system? 

We can think of initial data for Einstein's equation as given a picture of the 
space time at a given time. It consists of a three dimensional manifold S with 
some fields on it. The fields must satisfy the constraint equations. If the data 
describe an isolated system, the manifold S is naturally divided in two regions. 
One compact region fl which "contains the sources" , and its complement S\Q 
which is unbounded. We will call the latter one the asymptotic region. Of course, 
Q can be as large as we want, the only requirement is that it be bounded. The 
fact that S \ Q is unbounded means that we can go as far a we want from the 
source region Q : this capture the idea that the sources in fl are isolated. That 
"the sources are in fl" suggest that the field decays in S \ fl. Then, the initial 
data approach flat initial data near infinity. We call them asymptotically flat 
initial data. An opposite situation is when we want to study the universe as 
a whole. In this case one usually consider initial data where the manifold S is 
compact without boundary. 

The simplest and most important example of asymptotically flat initial data 
is the case when S is M 3 and fl some ball. Such data can describe, for example, 
ordinary stars. Consider an initial data set for a binary system, as is shown in Fig. 
[l]. The stars are inside fl. Ordinary stars emit light, then in the asymptotic region 
we will have electromagnetic field beside the gravitational one. We can also have 
some dilute gas in this region. It is important to recall that we do not require 
vacuum in the asymptotic region. We only require that the field decays properly 
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Fig. 1. Initial data for two stars. The region SI contains the matter sources, which are 
represented by black regions. The manifold S is R 3 . 



outside SI. However, in many situations it is a good approximation to assume 
vacuum in that region. When they evolve, the stars will follow some trajectory. 
Presumably an spiral orbit. The radius of the orbit will decrease with time, since 
the system loses energy in the form of gravitational waves. At late times, the 
system will settle down to a final stationary regime. This final state can be a 
rotating, stationary star. Or, if the initial data have enough mass concentrated 
in a small SI, a spinning black hole. We show the conformal diagrams of this 
two cases in Fig. and Fig. [| The gravitational radiation is measured at null 
infinity, where the observer is placed. 
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Fig. 2. Conformal diagram of the space 
time generated by the initial data given 
in Fig. [l], in the case in which the hnal 
state is a stationary star. 



Fig. 3. Conformal diagram of the space 
time generated by the initial data given 
in Fig. in the case in which the final 
state is a black hole. 



Since the sources are in SI, the fields in this region can be very complicated. 
They depend on the kind of matter that forms the stars. Remarkable, even the 
topology of SI can be complicated. The example in Fig. |l| has trivial topology, 
however non trivial topologies are relevant for black hole initial data. Consider 
the initial data for the Schwarzschild and Kerr black hole. In this case S = 5 2 xl 
and SI — S 2 x [—a, a] for some constant a, as it is shown in Fig. ||. The asymptotic 
region of these data has two disconnected components S \ SI = Si + S^. Each 
component is diffeomorphic to R 3 minus a compact ball. We say that the data 
have two asymptotic ends. That is, there exist two disconnected "infinities". 
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Observations can be made in one of them but not in both. It is not clear the 
physical relevance of this extra asymptotic end. A normal collapse, as is shown 
in Fig. I], will not have it. Therefore, one can think that only one asymptotic end 
is astrophysically relevant, being the other one just a mathematical peculiarity. 
The topology is not fixed by the Einstein equations. In this sense it remains quite 
arbitrary, one can thing that it is a kind of boundary conditions that has to be 
extra imposed to the equations. However, non trivial topologies appear naturally 
in the study of vacuum stationary black holes. A space time is stationary if it 
admits a timelike Killing vector field. The Schwarzschild and Kerr metrics are 
stationary. One can prove that every vacuum stationary asymptotically flat space 
time with trivial topology must be Minkowski. That is the non trivial topology 
is the "source" of the gravitational field in the vacuum stationary space times. 
In Fig. we show another example with different a topology, in this case we 
have only one asymptotic end. Other topologies with many asymptotic ends are 
of interest because their evolution in time may describe the collision of several 
black holes. In Fig. || and Fig. [ll] we show some examples. 



One usually has the idea that the matter sources generate gravitational field. 
To some extend this is true for many cases of astrophysical interest. However, 
as we have seen, a pure vacuum initial data with non trivial topologies can have 
mass. Moreover, one can also have a pure gravitational radiation initial data, 
that is a vacuum initial data with trivial topology with non zero mass. These 
kind of data are important in order to understand properties of the radiation 
itself which do not depend on the specific matter models. It is even possible to 
form a black hole with these type of data (see 0). 

In contrast to f2, the asymptotic region is very simple. The fields on it are 
approximately flat and its topology is M 3 minus a compact ball. We want to 
analyze the fields in this region. It is of course true that the fields there are 
determined by the fields in Q. But some important properties of them do not 
depend very much on the detailed structure of the sources in Q. The most 
important of these properties is that the initial data will have a positive mass if 
the sources satisfies some energy condition. In this article we want to study some 




Fig. 4. One black hole initial data. 
The asymptotic region has two disjoint 
components Si and S2- The compact 
set Q is represented by the shadowed 
region. 



Fig. 5. Initial data for a Misner worm- 
hole The asymptotic region has 
only one component. The compact set 
O is represented by the shadowed re- 
gion. 
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other properties of the fall of the fields in the asymptotic region which will be 
also independent on the sources in ft. It is important to note that we can only 
measure precisely this kind of properties, because we can not prepare the initial 
conditions of an astrophysical system in the laboratory and then it is impossible 
to find the exact initial data for a real system. We can only analyze those effects 
which do not depend very much on fine details of the sources. 

In order to simplify the hypothesis of the theorems, we will assume in this 
article that there are no matter sources in S. All the results will remain valid 
under suitable assumptions on the matters sources. 

We summarize the discussion above in the following definition. An initial data 
set for the Einstein vacuum equations is given by a triple (S, hab^ab) where S is 
a connected 3-dimensional manifold, h a b a (positive definite) Riemannian metric, 
and ^ab a symmetric tensor field on S. They satisfy the vacuum constraint 
equations 

b b § ab - b a 4> = o, (i) 

R + § 2 - ^ ab ^ ab = 0, (2) 

on S, where D a is the covariant derivative with respect to h a b, R is the trace 
of the corresponding Ricci tensor, and 'P = h^'Pab- The data will be called 
asymptotically flat with N asymptotic ends, if for some compact set Q we have 
that S\Q = Ylk=i ^ fe ' where §k are open sets such that each Sk can be mapped 
by a coordinate system i? diffeomorphically onto the complement of a closed 
ball in R 3 such that we have in these coordinates 

9m 

h t3 = (l + —)6 l3 +0(?- 2 ), (3) 

^=0(r- 2 ), (4) 

as f = Q3j=i(^"') 2 ) 1/ ' 2 ~~ * 00 m eac h set Sk- Here the constant m denotes the 
mass of the data, a, b, c... denote abstract indices, i,j,k..., which take values 
1,2,3, denote coordinates indices while 5ij denotes the flat metric with respect 
to the given coordinate system x J . Tensor indices will be moved with the metric 
h a b and its inverse h ab . We set Xi = x % and d l = <9;. These conditions guarantee 
that the mass, the momentum, and the angular momentum of the initial data 
set are well defined in every end. 

We want to analyze the higher order terms in (||) and (^). For example, the 
terms (3(f~ 2 ) could have the form f~ 2 sin(f 10 ). This function is certainly 0(f~ 2 ), 
but any derivative of it will blow up at infinity. Should such complicated terms 
be admissible in a description of a realistic isolated systems? In Newton's theory 
and Elcctromagnetism one can give a definitive answer to this question: the fields 
have a fall off behavior like powers in a radial coordinate. Take for example a 
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matter density p with compact support. The Newton's gravitational potential u 
satisfies the Poisson equation 

Au = 4irp, (5) 

where A is the Laplacian. Outside the support of p the potential u satisfies 
Au = 0. It is a well known that every harmonic function that goes to zero at 
infinity has an expansion of the form 



u k - x 



f.k ' 



(0) 



k>l 



where u k are the spherical harmonics of order k. In this case the field equations 
force the potential to have the fall off behavior (Q). 

The situation for the Einstein equation is more complicate. In analogy with 
(||) , one can ask the question whether there exist a class of initial data such that 
the metric and the extrinsic curvature admit near space-like infinity asymptotic 
expansions of the form 



f 



fy-a+^to+E^' (7) 



k>2 



ft 



E 

fe>2 



\f/ k 



(8) 



where h k j and are smooth function on the unit 2-sphere (thought as being 
pulled back to the spheres f — const, under the map x J — > x 3 /r). In this article I 
want to give an answer to this question. It is not only for convenience or aesthetic 
reasons that we want to avoid terms like f _2 sin(f 10 ) in the expansions. They 
are also very difficult to handle numerically. 

In order to see how the gravitational field behaves near infinity, it is natural 
to study first some examples. Consider asymptotically flat static space-times. 
We say that a space-time is static if it admits a hypersurface orthogonal time 
like Killing vector. One can take one of this hypersurfaces and analyze the fields 
h a i, and ft^ on it. The simplest static space time is the Schwarzschild metric. 
In this case we have hij = (1 + m/(2f)) 4 5y, ftj = 0, which is certainly of the 
form (0), (|8|). For general static, asymptotically flat space times it can be proved 
that the initial data have also asymptotic expansions of the form (0) and (|8|). 
Moreover, in this case the fields are analytic functions of the coordinates. Then 
the expansions (0) and (||) are in fact convergent powers series, in analogy with 

t. This result is far from being obvious. It was proved by Beig and Simon in 
based on an early work by Geroch . 

One can go a step further and consider asymptotically flat stationary solu- 
tions. A stationary space time admits a timclike Killing field which is in general 
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non hypersurface orthogonal. These space times describe rotating stars in equi- 
librium. The twist of the Killing vector is related to the angular velocity of the 
star. In this case there are no preferred hypersurfaces. The most important ex- 
ample of stationary solution is the Kerr metric. It can be seen that for t = const. 
slice of the Kerr metric in the standard Boyer-Lindquist coordinates the fields 
h a b and & a b satisfy (]?]) and (||). For general stationary asymptotically flat so- 
lutions there also exist slices such that (Q) and (|^) holds. The essential part of 
this result was also proved by Beig and Simon in (see also However, 
in these works the expansions are made in the abstract manifold of the orbits 
of the Killing vector field. In contrast with the static case this manifold does 
not correspond to any hypersurface of the space time. One has to translate this 
result in terms of the metric and the extrinsic curvature of some slices. This last 
step was made in [Q. As in the static case, the fields are analytic functions of 
the coordinates. 

We have shown that there exist non trivial examples of initial data which 
satisfies (0) and (0). But how general are these expansions? For example, is it 
possible to have data with gravitational radiation that satisfies (0) and (||)? I 
want to show that in fact there exists a large class of asymptotically flat initial 
data which have the asymptotic behavior (Q) and (^J). These data will not be, in 
general, stationary. 

The interest in such data is twofold. First, the evolution of asymptotically 
flat initial data near space-like and null infinity has been studied in considerable 
detail in [ p2[ . In that article has been derived in particular a certain "regularity 
condition" on the data near space-like infinity, which is expected to provide a 
criterion for the existence of a smooth asymptotic structure at null infinity. To 
simplify the lengthy calculations, the data considered in 22 have been assumed 
to be time-symmetric and to admit a smooth conformal compactification. With 
these assumptions the regularity condition is given by a surprisingly succinct 
expression. With the present work we want to provide data which will allow 
us to perform the analysis of without the assumption of time symmetry 
but which are still "simple" enough to simplify the work of generalizing the 
regularity condition to the case of non-trivial second fundamental form. Second, 
the "regular finite initial value problem near space-like infinity" , formulated and 
analyzed in j22], suggests how to calculate numerically entire asymptotically 
flat solutions to Einstein's vacuum field equations on finite grids. In the present 
article I provide data for such numerical calculations which should allow us to 
study interesting situations while keeping a certain simplicity in the handling of 
the initial data. 

The difficulty of constructing data with the asymptotic behavior (0), (3) 
arises from the fact that the fields need to satisfy the constraint equations (1) 
and Part of the data, the "free data", can be given such that they are 
compatible with ([?]), (^). However, the remaining data are governed by elliptic 
equations (the constraint equations will reduce to elliptic equations as we will 
see below) and we have to show that (Q), (||) are in fact a consequence of the 
equations and the way the free data have been prescribed. 
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To employ the standard techniques to provide solutions to the constraints, 
we assume 

9 = 0, (9) 

such that the data correspond to a hypersurface which is maximal in the solution 
space-time. 

I give an outline of the results which are available so far (see [lis] and ]l5| for 
more details and proofs). Because of the applications indicated above, we wish 
to control in detail the conformal structure of the data near space-like infinity. 
Therefore we shall analyze the data in terms of the conformal compactification 
(S, h a b, tf'ab) of the "physical" asymptotically flat data. Here S denotes a smooth, 
connected, orientable, compact 3-manifold. Take an arbitrary point i in S and 
define S = S\{i}. The point i will represent, in a sense described in detail 
below, space-like infinity for the physical initial data. The physical manifold S is 
essentially the stereographic projection of the compact manifold S. For example 
if we chose S to be 5* 3 then S will be R 3 . 




Fig. 6. The stereographic projection of S 3 in R 3 . A neighborhood B a of the point i is 
mapped into the asymptotic region of the physical manifold S. 



Working with S and not with 5, has several technical advantage. It is simpler 
to prove existence of solutions for an elliptic equation on a compact manifold 
than in a non compact one. The price that we have to pay is that the equations 
will be singular at i. However this singularity is mild. It is also simpler to analyze 
the fields in terms of local differentiability in a neighborhood of i than in terms 
of fall off expansions at infinity in S. Moreover, this technique is also useful to 
construct initial data with non-trivial topology. By singling out more points in 
S and by treating the fields near these points in the same way as near i we could 
construct data with several asymptotically flat ends. In Fig. |?] -|l3|we show some 
examples. All the following arguments equally apply to such situations, however, 
for convenience we restrict ourselves to the case of a single asymptotically flat 
end. 

We assume that h a {> is a positive definite metric on S with covariant derivative 
D a and 9 a b is a symmetric tensor field which is smooth on S. In agreement with 
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Fig. 9. The Brill-Lindquist initial data 
for two black holes . The data have 
three asymptotic ends. 




Fig. 11. The Misner initial data for 
two black holes p^ ]. The data have 
three asymptotic ends. 




Fig. 10. A compactification of these 
data in the sphere. 




Fig. 12. A compactification of these 
data in the torus. 



Fig. 13. The Misner wormhole initial Fig. 14. A compactification of these 

data. The data have only one asymp- data in the torus, 

totic end. 

(||) we shall assume that ^/ a b is trace free, 

h ab V ab = Q. 

The fields above are related to the physical fields by rescaling 

Kb = o A h ab , # ab = e- 2 V ab , (10) 

with a conformal factor 9 which is positive on S. For the physical fields to satisfy 
the vacuum constraints we need to assume that 

D a V ab = on S, (11) 
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(D b D b - -R)6 = --V ab V ab 0- 7 on S. (12) 
o o 

Equation ( [l2|) for the conformal factor 9 is the Lichnerowicz equation, transferred 
to our context. 

Let x 3 be /i-normal coordinates centered at i such that h^i — Ski at i and set 
r = {J^i^ii^) 2 ) 1 ^ 2 ■ To ensure asymptotical flatness of the data ( |Io|) we require 

<f r ab = 0(r^ 4 ) as r^O, (13) 



limr0 = l. (14) 



-o 



In the coordinates x j = x 3 /r 2 the fields © will then satisfy @, (@) (cf. |2J 
[ p2| for this procedure). 



Not all data as given by fllO|), which are derived from data h a b, ^ab as de- 
scribed above, will satisfy conditions (Q), (||). We will have to impose extra 
conditions and we want to keep these conditions as simple as possible. 

Since we assume the metric h ab to be smooth on S, it will only depend on 
the behavior of near i whether condition (fjj) will be satisfied. Via equation ( p"2| ) 
this behavior depends on & ab . What kind of condition do we have to impose on 
\P ab in order to achieve (0) ? 

The following space of functions will play an important role in our discussion. 
Denote by B a the open ball with center i and radius r = a > 0, where a is 
chosen small enough such that B a is a convex normal neighborhood of i. A 
function / 6 C°°(5') is said to be in E°°{B a ) if on B a we can write f — fi+ rf 2 
with /i,/a £ C°°(B a ) (cf. definition |l|). An answer to our question is given by 
following theorem: 

Theorem 1. Let h ab be a smooth metric on S with positive Ricci scalar R. 
Assume that W ab is smooth in S and satisfies on B a 

r 8 W ab & ab e E°°(B a ). (15) 

Then there exists on S a unique solution 6 of equation which is positive, 

satisfies ftlA), and has in B a the form 



9 = -, ieB°°(B a ), ${i) = l. (16) 
r 

In fact, we will get slightly more detailed information. We find that 6 = 
ui+ru2 on B a with u 2 & E°°(B a ) and a function u\ <E C°°(B a ) which satisfies 
ui = 1 + 0(r 2 ) and 

{D b D b - ^R)-y = Qr-, 
in B a \{i}, where 9r is in C°°(B a ) and vanishes at any order at i. 
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If 6 has the form ( |l6| ) then (^) will be satisfied due to our assumptions on 
h a b- 

Note the simplicity of condition (|l5|). If the metric is analytic on B a it can be 
arranged that 6^ =0 and u% is analytic on B a (and unique with this property, see 
[p3[). The requirement R > 0, which ensures the solvability of the Lichnerowicz 
equation, could be reformulated in terms of a condition on the Yamabe number 
(cf. f§).' 

Theorem |l| has two parts. One is the existence and uniqueness of the solution. 
This part depends on global properties of the fields on S. It can be proved under 
much weaker assumptions on the differentiability of h a b and <P a b. The second part 
concerns the regularly in B a , and depends only on local properties of the fields in 
B a . This part can also be proved under weaker hypothesis. These generalizations 
have physical relevance, I will come back to this point in the final section. 

It remains to be shown that condition (|l^) can be satisfied by tensor fields ^ a b 
which satisfy (pd|), (|l3|). A special class of such solutions, namely those which 
extend smoothly to all of S, can easily be obtained by known techniques (cf. 
[p~4| ) . However, in that case the initial data will have vanishing momentum and 
angular momentum. To obtain data without this restriction, we have to consider 
fields if'ah e C°°(S) which are singular at i in the sense that they admit, in 
accordance with (El), (lOh, (M), at i — {r — 0} asymptotic expansions of the 



form 



~ ®if k with !pg G C°°(S' 2 ). (17) 

fc>-4 

It turns out that condition ( |l5| ) excludes data with non- vanishing linear mo- 
mentum, which requires a non- vanishing leading order term of the form 0(r~ 4 ). 
In section ^| we will show that such terms imply terms of the form log r in 9 and 
thus do not admit expansion of the form (^). However, this does not necessarily 
indicate that condition ([l5]) is overly restrictive. In the case where the metric h a b 
is smooth it will be shown in section [2] that a non- vanishing linear momentum 
always comes with logarithmic terms, irrespective of whether condition ( p"5| ) is 
imposed or not. 

There remains the question whether there exist fields ^ab which satisfy ( p"5| ) 
and have non-trivial angular momentum. The latter requires a term of the form 
(3(r~ 3 ) in (|l7j). It turns out that condition ( |l5|) fixes this term to be of the form 

= 4r (3n a n b - 5 ab ) + ^(n a e bcd J c n d + n b e acd J c n d ), (18) 

where n l = x l jr is the radial unit normal vector field near i and J k , A are 
constants, the three constants J k specifying the angular momentum of the data. 
The spherically symmetric tensor which appears here with the factor A agrees 
with the extrinsic curvature for a maximal (non-time symmetric) slice in the 
Schwarzschild solution (see for example ||). Note that the tensor W^j satisfies 
condition ([l5]) and the equation d a W^j = on S for the flat metric, hence it is 
a non-trivial example. We want to study more general situations. 
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In the existence proof for equations ([ll]) and (|l^), the possible existence of 
conformal Killing vectors £ a of the metric h a i will play an important role. A 
conformal Killing vector is a solution on S of £h£, = 0, where we have defined 
the conformal Killing operator 

{C h v) ab = D a v b + D b v a --h ab D c v c . (19) 

3 

Given £ a we define the followings constants 

S a = e\ c D b eW, a=±D a £ a (i). (20) 

The free-data in the solution <P ab consist in two pieces: a singular and a 
regular one, which we will denote by &j$ ng and <P^ b g respectively. We define 
<P a J ng in B a \ W by 

^ na = X^fj-\h cd Wfjh ab ), (21) 

and vanishes elsewhere. Here \ denotes a smooth function of compact support 
in B a equal to 1 on B a /2- We assume that <?°eg can be written near i in the form 

< s = r- 3 < es + < e9 , (22) 

where , ^reg are smooth symmetric trace free tensors in S and such that 
$Y reg = 0(r 2 ), and x t x } <P\ j reg = r 2 <P with some <PeC°°(B a ). 

Theorem 2. Let h a b a smooth metric in S. Let $> ab ing and @reg given- by ([f 
and /pT 



i) If the metric h admits no conformal Killing fields on S, then there exists 
a unique vector field v a such that the tensor field 

* ab = $% ng + $? eg + (£hvT b , (23) 

satisfies the equation D a 'L rab = in S\ {i}. The solution satisfies fy a b — 0(r~ 3 ) 
at i and condition nil 



ii) If the metric h admits conformal Killing fields £ a on S, a vector field v a as 
specified above exists if and only if the constants J a , A (partly) characterizing 
the tensor field &gi ng satisfy the equation 

J a S a +Aa = 0, (24) 

for any conformal Killing field £ a of h, where the constants S a and a are given 
by M). 



In both cases the angular momentum of <J rab is given by J a . This quantity 
can thus be prescribed freely in case (i). 
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Theorem [^, as Theorem |J contains two parts. One is the existence and 
uniqueness of the solution v a . The restriction (24) appears in this part. See 
[El, jl8| and |l^] for a generalization and interpretation of this condition. The 
other part concerns the regularity of the solution in B ai namely that it satisfies 
condition Jl5|). This part depends only on local assumptions of the fields in B a . 



We prove a more detailed version of this theorem in 1 18 



2 Solution of the Hamiltonian Constraint with 
Logarithmic terms 

Assume, for simplicity, that the metric h a b is conformally flat in B a . Then the 
operator that appears in the left hand side of equation (|l^) reduce to the Lapla- 
cian A. This simplification is of course a restriction on the allowed initial data, 
but it already contains all the main problems and the essential features of the 
more general case. It is also important to remark that even under this assump- 
tion it is possible to describe a rich family of initial data, since we are making 
restrictions only on B a . Consider the Poisson equation (^J). The example 

A{\ogrh m ) = r- 2 h rn (2m + 1), (25) 

where h m is an harmonic polynomial of order m (i.e. Ah m — 0) shows that 
logarithmic terms can occur in solutions to the Poisson equation even if the 
source has only terms of the form r s p where p is some polynomial. That is, even 
if our free data have expansion in powers of r logarithmic terms will appear in the 
conformal factor 9. We shall use this to show that logarithmic terms can occur in 
the solution to the Lichnerowicz equation if the condition r 8 ^ r a b^ rab £ E°°(B a ) is 
not satisfied. Our example will be concerned with initial data with non- vanishing 
linear momentum. 

We assume that in a small ball B a centered at i the tensor 'P ab has the form 

yab = gab + W ab^ ( 26 ) 

where is given in normal coordinates by 

= JL (-P V - P b n a - (6 ab - bn a n b ) P c n c ) , (27) 

and iP'jj' = 0(r~ 3 ) is a tensor field such that 1 & ab satisfies equation ( pTTj) . We will 
assume also that tf^ b satisfies some mild smoothness condition (cf. jjfq]). 

Since >Fp b is trace-free and divergence-free with respect to the flat metric, 
we could, of course, choose h a b to be the flat metric and l P^ b = 0. This would 
provide one of the conformally flat initial data sets discussed in |ic|| . We are 
interested in a more general situation. 

Lemma 1. Let h a b be a smooth metric, and let ^ ab be given by $2(\). Then, 
there exists a unique, positive, solution to the Hamiltonian constraint (f?j[). In 
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B a it has the form 

= ^ + \m + r ((9 (P'n;) 2 - 33 P 2 ) 

+ T6 mr2 (iP' + liHP'nif-P^logr^+un, (28) 

where P l Pi = P 2 , the constant m is the total mass of the initial data, w\ is a 
smooth function with wi = 1 + 0(r 2 ), and ur € C 2 ' a {B a ) with ur(Q) = 0. 

Since w\ is smooth and ur is in C 2a (B a ), there cannot occur a cancellation 
of logarithmic terms. For non-trivial data, for which m ^ 0, the logarithmic term 
will always appear. In the case where h a b is flat (in this case we have w± = 1) 
and \Pr = an expansion similar to ( |28| ) has been calculated in |^6| . 

3 Explicit solutions of the momentum constraint 

Instead of given the proof of Theorem EL in this section I want to present some 
explicit solutions of equation ( pd| ) for conformally flat and also for axially sym- 
metric metrics. In the first case we will construct all the solutions, in the second 
one only some of them. We will show how to achieve condition ( |l5| ) in terms of 
the free data. In both cases the solution is constructed in terms of derivatives 
of some free functions. That makes them suitable for explicit computations. In 
particular, if we assume that these free functions have compact support, then 
<P ab will also have compact support. We note incidentally that, as an applica- 
tion, one can easily construct regular hyperboloidal initial data with non-trivial 
extrinsic curvature. 

3.1 The momentum constraint on Euclidean space 

In the following we shall give an explicit constructing of the smooth solutions 
to the equation d a ^ ab = on the 3-dimensional Euclidean space E 3 (in suitable 
coordinates M 3 endowed with the flat standard metric) or open subsets of it. 
Another method to obtain such solutions has been described in Q, multipole 
expansions of such tensors have been studied in || . 

Let i be a point of E 3 and x k a Cartesian coordinate system with origin i 
such that in these coordinates the metric of E 3 , denoted by S a b, is given by the 
standard form Ski- We denote by n a the vector field on E 3 \ {i} which is given 
in these coordinates by x k /\x\. 

Denote by m a and its complex conjugate fh a complex vector fields, defined 
on E 3 outside a lower dimensional subset and independent of r = \x\, such that 

m a m a = fh a rh a = n a m a — n a rh a = 0, m a fh a = 1. (29) 

There remains the freedom to perform rotations m a — > e l *m a with functions / 
which are independent of r. 
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The metric has the expansion 

Sab = n a,nb + m a m b + m a fh bl 
while an arbitrary symmetric, trace-free tensor W ab can be expended in the form 

r 3 ^ab = £,(3n a n b - S ab ) + V2?7in (a m b) + V2fjin {a m b) + 

p,2m a mb + ^2m a fh b , (30) 

with 

£ = \r z V ab n a n\ m = V2r 3 W ab n a m b 1 ^ = r 3 V ab m a m b . 

Since W ab is real, the function £ is real while 771, fj.2 are complex functions of spin 
weight 1 and 2 respectively. 
Using in the equation 

d a ^ ab = 0, (31) 

the expansion ( |30"1) and contracting suitably with n a and m a , we obtain the 
following representation of (||l]) 

4rd r £, + Brji + 377! = 0, (32) 

r<9 r 77i + 3/j 2 - 3£ = 0. (33) 

Here d r denotes the radial derivative and 3 the edth operator of the unit two- 
sphere (cf. [j33| for definition and properties). By our assumptions the differential 
operator 3 commutes with d r . 

Let s Yi rn denote the spin weighted spherical harmonics, which coincide with 
the standard spherical harmonics Yi m for s — 0. The s Y; m are eigenfunctions of 
the operator 33 for each spin weight s. More generally, we have 

wrw.) =(-d- (i ( ';: ) - ) ," ( ;;t;" .^. <*> 

If fi s denotes a smooth function on the two-sphere of integral spin weight 
s, there exists a function /i of spin weight zero such that 77^ = 3 s 77. We set 
7] R = Re(?7) and ij 1 — i Im(?7), such that 77 = 77 + I) 1 , and define 

77f =3V, % 7 = 3V, 
such that t] s = 77^ + 77^ • We have 



3 S ?7« = d s r) R , 3V = -3 S 77 J 



Using these decompositions now for 771 and /J2, we obtain equation (32) in 
the form 



2rd r £ = — 33?7 



(35) 
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Applying 8 to both sides of equation ( |33| ) and decomposing into real and imag- 
inary part yields 

rd^n 1 = -8 2 SV (36) 
2rd r (rd r £) + 89? = 9 2 3V fl - (37) 



Since the right hand side of (35) has an expansion in spherical harmonics 
with I > 1 and the right hand sides of ([36]), (|37|) have expansions with I > 2, we 
can determine the expansion coefficients of the unknowns for I = 0, 1. They can 
be given in the form 

£ = A + rQ + -P, rj 1 = iJ + const., ij R = rQ - -P + const., 

r ' r 

with 

P='^P a n a , Q=^Q a n a , J = 3J a n a , (38) 



where A,P a ,Q a ,J a are arbitrary constants. Using (p0|), we obtain the corre 
sponding tensors in the form (cf. (|lC|]) 

3 
2~7 



Wf = {-P a n b - P b n a - (6 ab - 5n a n b ) P c n c ) , (39) 



*f = ^(n a e bcd J c n d + n b e acd J c n d ), (40) 

We assume now that £ and rj 1 have expansions in terms of in spherical har- 
monics with I > 2. Then there exists a smooth function A2 of spin weight 2 such 
that 

Using these expressions in equations ( |35"| ) - ( |37| ) and observing that for smooth 
spin weighted functions u s with s > we can have 3^t s = only if fj, s = 0, we 
obtain 

Bn R = -2rd r BX R , d 2 u' = -r<9 r A^, 

d 2 fi R = 2rd r (rd r X R ) - 2\ R + m\ R . 

We are thus in a position to describe the general form of the coefficients in the 
expression ( f30| ) 

£ = W\ R + A + rQ + ^-P, (43) 

r?i = -2rd r B\ R + d\i + rdQ- -dP + iBJ, (44) 

r 

li 2 = 2rd r {rd r \ R )-2\ R + m\ R -rdrXl (45) 
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Theorem 3. Let A be an arbitrary complex C°° function in B a \ {i} C E 3 with 
< a < oo, and set A 2 = 3 2 A. Then the tensor field 

yab = yab + yah + yah + yab + yab ^ ^gj 

where the first four terms on the right hand side are given by ( p?^ ) - while 
<Z^ b is is obtained by using in ( f%[ ) only the part of the coefficients fyty - 
which depends on A2, satisfies the equation D a W a b = in B a \ Conversely, 
any smooth solution in B a \ {i} of this equation is of the form 

Obviously, the smoothness requirement on A can be relaxed since tf^ b S 
C 1 (B a \ {i}) if A S C 5 (B a \ {i}). Notice, that no fall-off behavior has been 
imposed on A at i and that it can show all kinds of bad behavior as r — > 0. 

Since we are free to choose the radius a, we also obtain an expression for the 
general smooth solution on E 3 \ {i}. By suitable choices of A we can construct 
solutions if - ^ 6 which are smooth on E 3 or which are smooth with compact sup- 



port. Finally we provide tensor fields which satisfy condition (15) and thus prove 
a special case of theorem [| see [|l8| for the proof. 



Theorem 4. Denote byW ab a tensor field of the type (f^j). IfrX e E°°(B a ) and 
P a = 0, then r % V ab V ab £ E°°(B a ). 

We wish to point out a further application of the results above. Given a 
subset S of K 3 which is compact with boundary, we can use the representation 
(E™ to construct hyperboloidal initial data ((p0[) on S with a metric h which is 



Euclidean on all of S or on a subset U of S. In the latter case we would require !?"° 
to vanish on S\U . In the case where the trace-free part of the second fundamental 
form implied by h on dS vanishes and the support of & ab has empty intersection 
with dS the smoothness of the corresponding hyperboloidal initial data near the 
boundary follows from the discussion in (Q). Appropriate requirements on h 
and <P ab near dS which ensure the smoothness of the hyperboloidal data under 
more general assumptions can be found in (El). 



3.2 Axially symmetric initial data 

The momentum constraint with axial symmetry has been studied in H , |ll[ and 
[ [TGI . Assume that the metric h ab has a Killing vector rf , which is hypersurface 
orthogonal. We define 77 by 77 = rfr^h^. Following ^7j, consider 'P ab defined by 

V 

where ^ a satisfies 

£ v & a = Q, & a r] a = Q, D a V a = 0, (48) 
with the Lie derivative £ n with respect rf". 
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We use the Killing equation Di a r)b) = 0, the fact that ij a is hypersurface 
orthogonal, (i.e.; it satisfies D a r]b = —r}[ a D b ] In 77) and equations ( p8| ) to conclude 
that ^ ah is trace free and divergence free with respect to the metric h a b- 

The solution of equations (E8h can be written in terms of a scalar potential 



77 

Using this equation we find that 

1 Dm D c uj 
W ab W ab = 2^- 2 . 50 

We want to find now which conditions we have to impose in uj in order to achieve 
([Isj). The metric h a b has the form 

h, a b = e a b-\ , (51) 

V 

where e ab is the two dimensional metric induced on the hypersurfaces orthogo- 
nal to r) a , its satisfies £ ri e a i 3 = 0. All two dimensional metrics are locally con- 
formally flat, we will assume here that e ab is globally conformally flat. Then, 
we can perform a conformal rescaling of h a b such that in the rescaled metric 
the corresponding intrinsic metric e a b is flat. We will denote this rescaled metric 
again by h ab . Assume that rj a is a rotation. Take spherical coordinates (r, j5, 4>) 
such that rf — (d/d(f) a . In this coordinates the metric has the form 

h = dr 2 +r 2 di9 2 +r)d(/) 2 . (52) 

The norm 77 can be written as 

77 = r 2 sin 2 i9(l + /), (53) 

where / satisfies 1 + / > 0. Note that r is the geodesic distance with respect to 
the origin. 

In these coordinates equation (^) has the form 

2 (rW + (W) 
(1 + /) r° sm iSI 

In the flat case (i.e. / = 0) this solution reduce to the one given by Theorem 
I with 

P = Q = A = 0, Af = 0, (55) 

and 



iuj = \' 2 sin 2 i? + iJ 2 (-3cos79 + cos 3 '!?), 



(56) 
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where J z is the only non vanished component of the vector J a , and £^\\ = 0. 
Motivated by this expression we write to in the form 

lu = 3 2 Asin 2 7? + J z (-3cosi9 + cos 3 tf), (57) 

where A is an arbitrary zero spin real function which depends on r and i?. The 
constant J z will give the angular momentum of the data. This can be seen from 
the following expression 

JZ = ^ J s <Pan a dS, (58) 

where S is any closed two-surface in the asymptotic region and n a its normal. 
Using equation ( |57| ) and lemma 9 of fr§[ ] we prove the following lemma. 

Theorem 5. Assume that f g E°°(B a ). Let<P ab be given by @, @ and @, 
where J z is an arbitrary constant. IfrX G E°°(B a ) then r % V ab $ ab g E°°{B a ). 

Note that in Theorem || we have assumed only / g E° a (B a ), that is /, and 
hence the conformal metric h a b, is not required to be smooth. I will come back 
to this point in the final section. 



4 Main ideas in the proof of theorem |l| 

I want to describe in this section the main idea in the proof of theorem [|. I will 
concentrate on the regularity part of this theorem and not on the existence part, 
since the later is more or less standard. I will give an almost self contained proof 
of a simplified version of theorem |l|, which contains all the essential elements of 
the general proof. 

Consider the semi-linear equation 

Au = f(u,x), (59) 

on B a , where the function / is given by 

fM = - f {x) V(x)EE°°. (60) 
(1 + ru) ' 

This equation is similar to equation (|l2| ) when the metric h a b is flat in B a . 
Assume that we have a positive solution u g C 2 ' a (B a ). There exist several 
method to prove existence of solutions for semi-linear equations, see for example 
| fu)| for an elementary introduction to the subject. We want to prove the following 
theorem. 

Theorem 6. If u g C 2 ' a (B a ) and u > is a solution of equation (p9|), then 
ueE x (B a ). 
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The important feature of equation (J5^) is that the radial coordinate r ap- 
pears explicitly in /. As a function of the Cartesian coordinates x l , r is only in 
C a . Thus, we can not use standard elliptic estimates in order to improve the 
regularity of our solution. Instead of this we use the following spaces. 

Definition 1. For m G No and < a < 1, we define the space E m ' a (B a ) as 
the set E m ^ a {B a ) = {f = h + rf 2 : f u f 2 G C m > a (B a )}. Furthermore we set 
E°°(B a ) = {/ = h + rh ■ h,h G C°°(B a )}. 

The spaces E m,a has two important properties. The first one is given by the 
following lemma which is an easy consequence of definition [j]. 

Lemma 2. For f,g£ E m,a (B a ) we have 

(i) f + g G E m ' a (B a ) 

(ii) fg£E m > a (B a ) 

(in) Iff ^0 in B ai then l/f G E m > a (B a ). 
Analogous results hold for functions in E°°(B a ). 

The second important property is related to elliptic operators. Let u be a 
solution of the Poisson equation 

An = f. (61) 

Then we have the following lemma. 

Lemma 3. / G E m ' a (B a ) u G E m+2 > a (B a ). 

Proof: Since fa G C m,a (B a ) we can define the corresponding Taylor poly- 
nomial T m of order m. Define f 2 by f 2 — T m + f 2 . It can be seen that 
f2 = 0(r ra+a ). By explicit calculation we can prove that there exist a poly- 
nomial p m of order m such that A(r 3 p m ) = rT m . Set u = r p m + ur. Then ur 
satisfies the equation 

Au R = f 1 +rf 2 R . (62) 

One can prove that rf 2 G C rn ' a (B a ). This is not trivial because r is only in 
C a {B a ). Here we use that f£ = 0{r m+a ) (see lemma 3.6 in @). Then the right- 
hand side of equation ( |62| ) is in C m,a (B a ). By the standard Schauder elliptic reg- 
ularity (see H) we conclude that u fl G C m+2 < a {B a ). Thus u G E m+2 ' a (B a ). U 

As an application of lemmas and |^ we can prove theorem [| Using lemma 
||, we have that the function / satisfies satisfies the following property 

u G E m > a (B a ) =► /(n^), i) G i? m ' tt (Ba), (63) 

for every positive it. Assume that we have a solution u G C 2 ' a of equation (|59|). 
Using (|6^), lemma H and induction in m we conclude that u G E°°(B a ). Here we 
have used that if / G E m ' a (B a ) for all m G N 0) then / G E°°(B a ), see lemma 
3.8 in @. 

Note that the only property of / that we have used in order to prove theorem 
[| is ( |63| ) . Lemma [| can be generalized for second order elliptic operators with 
smooth coefficients (cf. [p"8|). 
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5 Final Comments 

In this section I want to make some remarks concerning the differentiability of 
the initial data. It is physically reasonable to assume that the fields are smooth 
in the asymptotic region S \ fl. However in the source region Q smoothness is 
a restriction. For example, the matter density of an star is discontinuous at the 
boundary. Generalizations of the existence theorems to include such situations 
have been made in |b| and (T^] . 

It is also a restriction to assume that the conformal metric is smooth in a 
neighborhood B a of the point i. The point i is the infinity of the data, hence 
there is no reason a priori to assume that the fields there will have the same 
smoothness as in the interior. Moreover it has been proved in ]l7j that the 
stationary space times do not have a smooth conformal metric. In this case the 
metric h a b in B a has the following form 

l>u ■ r :i lr r (64) 

where hjj and tif^ are analytic functions of the coordinates x 1 . Theorems |l| and || 
are generalized in jl^] for metrics that satisfy (Q). An example of this general- 
ization is Theorem g in which we have assumed that / 6 E°°, this assumption 
allows metrics of the form d&jj). 
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